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Abstract 

We study the 0(n) loop model on the honeycomb lattice with open boundary con- 
ditions. Reflection matrices for the underlying Izergin-Korepin .R-matrix lead to three 
inequivalent sets of integrable boundary weights. One set, which has previously been 
considered, gives rise to the ordinary surface transition. The other two sets correspond 
respectively to the special surface transition and the mixed ordinary-special transition. 
We analyse the Bethe ansatz equations derived for these integrable cases and obtain the 
surface energies together with the central charges and scaling dimensions characterizing 
the corresponding phase transitions. 
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1 Introduction 



The 0(n) model Q in two dimensions has attracted considerable attention in recent 
years (see, e.g. Refs. [2-21] and references therein), not least because of its connection to 
the problem of planar self-avoiding walks or polymers in the n — > limit |22| . Methods 
used to study the model have included, among others, Coulomb gas techniques || [|, ||, 
conformal invariance || |6|, [7|, [13], [15], |21|, Bethe ansatz [|, [It], [11], |T^, |19|| , bootstrap S- 
matrix [fEBl O, 181 EU1 and numerical transfer matrix calculations flifl . The self-avoiding 



walk problem has also been extensively studied via a number of techniques (see, e.g. [P3l 
and references therein). Of particular relevance here are the finite-size transfer matrix 
calculations |24], |25| . 

The critical behaviour of the 0(n) model is well understood, with an exact Bethe 
ansatz solution on the periodic honeycomb lattice ||, 1C , 1 1 providing a confirmation 
of the earlier Coulomb gas and conformal invariance based results in the bulk. However, 
at least from the perspective of exactly solved lattice models, the situation is not as 
satisfactory for surface critical behaviour. A Bethe ansatz solution on the open honey- 
comb lattice was found [[TJJ, again confirming earlier results R [], [23] for the ordinary 
surface transition P^| . However, predictions also exist [L^, [2Cj, |21| for other surface 
transitions [^6[. In this paper we derive exact results for the special surface transition, 
also on the honeycomb lattice, using a Bethe ansatz technique which utilises reflec- 
tion matrices in an essential way [p7| . We also derive exact results for the model with 
mixed ordinary and special boundary conditions. These results for the surface critical 
behaviour both complement and, in some cases, extend those currently available in the 
literature from other means. Our results for the special transition and its relevance to 



the problem of the adsorption transition for polymers were announced in |28[ . 

The 0(n) or n- vector model on the honeycomb lattice is a statistical mechanical 
model of n-dimensional unit vectors Si living on the sites i of the lattice and interacting 
via nearest neighbour couplings. Its partition function is given by 



Z' 



n d ^ 



exp 



(ij) 



On an infinite lattice the coupling constant is usually chosen to assume the same 
value J for all nearest neighbour pairs (ij). In this paper we study the model with 
open boundary conditions (see Fig 1) with taking values Ji (respectively, J r ) if the 
nearest neighbour pair (ij) lives on the left (right) boundary. Following Domany et. 
al. 0, we argue that its critical behaviour is in common with that of a more amenable 
(but still O(n)-symmetric) model whose partition function is given by 



z = [ (n d ^)n( i +^^-4- 



;i-2) 



i 



Figure 1: Honeycomb lattice of "width" iV = 8 on which the 0(n) model is 
defined. The lattice extends to infinity in both vertical directions. A typical 
configuration Q contributing to the partition sum is shown in bold, with 
P = 2, N t = 2, JV r = 4 and N b = 36. 

Here the nearest neighbour coupling constant Kij takes values Ki (respectively, « r ) if 
the pair (ij) lives on the left (right) boundary, and k otherwise. 

Through a "high-temperature expansion" (in small k, K\ and K r ) Z can be equiva- 
lently written as the partition sum 

Z = 5> J VV^n p , (1.3) 
Q 

over all configurations Q of P non-intersecting and closed loops of fugacity n on the 
same lattice, with N (respectively, N r ) bonds of Q living on the left (right) boundary 
and iV?, bonds not living on either boundary. When n — > only the empty graph 
contributes to the partition function. However, the i lines joining the two points x and 
fin a "watermelon correlator" G((x — y) survives, giving rise in the continuum limit 
to the correlation function of i self- and mutually avoiding walks tied at the ends. 

The partition sum (|1.3|) can also be written in the following way which will eventually 
allow a mapping of the (non-local) loop model onto a (local) vertex model: 

Z~£<i--.^n p (1.4) 

Q 

where Wi are Boltzmann weights for the allowed vertices depicted in Figure 2, of which 
there occurs Ni such in the configuration Q. The weights Wi are given by W\ — ■ ■ ■ — 

W 6 — 1, Wj — W 8 — = t, Wg — Wio = 1, It'll = KK^ 2 = U and W\2 = kk~ 2 = t r . The 

partition sum (|1.4| ) is defined up to an unimportant overall multiplicative factor. 
For the 0(n) model on the (infinite) honeycomb lattice, the critical line 

K- 1 = (2 + (2-n) 1 / 2 ) 1/2 (1.5) 
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Figure 2: Allowed loop vertices for the honeycomb loop model. Vertex i 
carries a Boltzmann weight tOj. 

was first obtained by Nienhuis || using a mapping to various models and the Coulomb 
gas. Baxter |J subsequently showed that a vertex model to which the loop model ( |1.3| ) 
can be mapped is solvable by the coordinate Bethe ansatz precisely along the line (|1.5|) . 
The Bethe ansatz solution was then used to derive critical properties for the 0{n) model 



in Refs. [10|, [Tl| • For instance, in the parametrization 



n = — 2cos7rg, (1.6) 

K ~ 1 = t = 2sin^, (1.7) 

the model has central charge 

6(1 — t;) 2 , . 

c=l- — — . 1.8 

9 



For the open honeycomb lattice of relevance to this paper, Batchelor and Suzuki [IS 
showed that along the line specified by Q1.5Q and, in addition, the model 

is likewise solvable by the coordinate Bethe ansatz and corresponds to the 0(n) model 
at the ordinary surface transition. The Bethe ansatz analysis showed that the central 
charge ( [L.8|) remains valid. Furthermore, the model has the set of scaling dimensions 



X°~° = h e+ltU (1.9) 
in agreement with an earlier identification in terms of the Kac formula, where |j] 

h m = W - x m + g2 ~ ( f ~ 1)2 - (l-io) 



Physically, this set of geometric dimensions characterizes the algebraic decay of the 
half-watermelon correlators Gf (x — y) at criticality H. 
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The present authors have recently presented a systematic approach [^7] to obtain 
integrable loop models with open boundaries. This approach begins with Sklyanin's 
extension of the Quantum Inverse Scattering Method to obtain integrable open spin 
chains p9| , pPf] , By incorporating inhomogeneities and making them alternate with the 
spectral parameter the Sklyanin transfer matrix becomes relevant to vertex models 
with open boundaries. For certain choices of .fT-matrices, or solutions to the reflection 
equation, as input to the Sklyanin transfer matrix the corresponding vertex model can 
be mapped to a loop model. 

The 0(n) model of interest here can be obtained from a special case of the Izergin- 
Korepin or A 2 vertex model |}3, [34| (see also the reference to Korepin in Ref. 
H). There are three inequivalent diagonal i^-matrices associated to the -R-matrix 
Of these the two giving rise to non-quantum group invariant spin chains allow a 



mapping to loop models. We showed in p7j how to recover the model of [19fl which 
corresponds to the ordinary surface transition. In addition we also found a different set 
of boundary weights for which the model is integrable; these correspond to 



-1/2 



Kl = K r = {-ZCOS^f) = K Sp . 

Note that this critical coupling is related to n through 

K Sp ={2-n)- 1 ' 4 



;i.i2) 



We shall presently argue that this solution corresponds to the 0(n) model at the special 
transition, with central charge ( |1.8| ) and the set of scaling dimensions 

^Sp-Sp 



x: 



h 



£+1,3, 



;i.i3) 



confirming a conjecture using bootstrap boundary S- matrices We also obtain a 

third set of integrable boundary weights which is a mixture of the first two; namely, 



K ■ ; ■ 



KSp- (1-14) 

We argue that this corresponds to the 0(n) model with mixed ordinary-special bound- 
ary conditions, with central charge 

6(2 -gf 



c 



O-Sp 



1 
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;i.i5) 



or equivalently c° _Sp = c — 24 t° _Sp with the "mixed-boundary" scaling index given 
by t°~ Sp = hi } 2, confirming a recent result obtained using conformal invariance [2l|. In 
addition, we find the set of scaling dimensions 



v O-Sp 



he+i,i — 2^ 



;i.i6) 



A brief summary of our results for mixed boundaries and their relevance to two- 



dimensional polymers has been given elsewhere [36 . 
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(2) 

2 Honeycomb limit of the A 2 loop model 



In this section we present the derivation of integrable 0(n) loop models on the hon- 

(2) 

eycomb lattice with open boundaries from the A 2 -R-matrix and associated diagonal 
i^-matrices |35|| . This involves the following steps: (I) Interpreting the Sklyanin trans- 
fer matrix t(u, u>) with alternating inhomogeneities JHj] w as a diagonal-to-diagonal 
transfer matrix to(^) for a vertex model on the square lattice with open boundaries. 
(II) Mapping the vertex model to a loop model on the same lattice. (Ill) Taking the 
'the honeycomb limit' to arrive at the honeycomb lattice loop model and its transfer 
matrix Via the same chain of mappings the eigenvalue expression for t(u, u>) (and 
associated Bethe ansatz equations) gives rise to the corresponding expression for £h- 
This procedure has been explained in detail in our previous paper |27j for several three- 
state i?-matrices. In particular, the 0(n) models with non-mixed boundaries to which 
we hitherto refer as 0-0 and Sp-Sp are treated therein. We therefore concentrate here 
on the mixed boundary case denoted by O-Sp. 



2.1 Sklyanin transfer matrix and eigenvalue expression 

The i?-matrix for the A 2 model |32[ is given in terms of a spectral parameter u and 
anisotropy A by 

\ 
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(2.1) 



where 



a{u) = sinh(u — 3A) — sinh(5A) + sinh(3A) + sinh(A), 

b{u) = sinh(w — 3A) + sinh(3A), 

c{u) = sinh(u — 5A) + sinh(A), 

d(u) = sinh(-u — A) + sinh(A), 

e{u) = -2e"t sinh(2A) cosh(f - 3A), 

ii 

e{u) = -2e2 sinh(2A)cosh(f -3A), 

f(u) = -2e- u+2A sinh(A) sinh(2A) - e~ A sinh(4A), 
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f(u) = 2e u ~ 2X sinh(A) sinh(2A) - e A sinh(4A), 

g(u) = 2e~2+ 2A sinh(f)sinh(2A), 

g(u) = -2e2- 2A sinh(|)sinh(2A). 

It was shown in that there exists three inequivalent diagonal .fT-matrices K~(u) 



satisfying the reflection equation 

1 2 
R 12 (u - v) K~ {u)R 2l {u + v) K~ (v) = 

2 1 

K~ (v)R 12 (u + v) K~ (u)R 21 (u - v). (2.2) 
They are namely, K~(u) = 1 and K~(u) = T ± (u), where 

^(u) = diag(e-> ± (u), ^(u), e^iv)), (2.3) 

with 

^±(u) = cosh(f -3A) ±isinh(|) (2.4) 

= cosh(f + 3A) =Fisinh(|). (2.5) 

Let K + (u) = K~(-u - 7]YM, where K~ ~(u) is any of the three if-matrices defined 
above, rj = — 6A — in is the crossing parameter and M = diag(e 2A , 1, e _2A ) is the crossing 



matrix. Then the Sklyanin transfer matrix [29, 30 



t(u, u) = tr a K + {u)T a (u, u) K~ (u)f a {u, w), (2.6) 

forms a commuting family: [t(u,u)),t(v,u))] = for arbitrary u, v. Here T a and T a are 
monodromy matrices defined by 

T a (u,u) = R al (u + ui)R a2 {u + uj 2 ) ■ ■ ■ RaN{u + uj n ), (2.7) 
f a (u,u)) = R Na (u-uj N )---R 2a (u-uj 2 )R la (u-uji), (2.8) 

where a labels the auxiliary space V a and 1, . . . , N label quantum spaces. With three 
choices each for K~{u) and K + (u) we have therefore nine possiblities for t{u,uS). The 
three non- mixed cases were dealt with in ]27|. It was shown there that the non quantum 



group-invariant cases, 0-0 and Sp-Sp, based on T ± (m) give rise to loop models, in 
contrast to the quantum group-invariant case with K~(u) = 1, K + (u) = M. Similar 
reasoning shows that only two of the mixed cases O-Sp and Sp-0 give rise to loop 
models. These two loop models are equivalent, as to be expected, and we henceforth 
deal only with the case O-Sp specified by K~(u) = T~(u) and K + {u) = T + (—it — rjY M . 
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The eigenvalue expression and Bethe ansatz equations for the Sklyanin transfer 
matrix t(u, uj) were obtained for the cases 0-0 and Sp-Sp in [27] using what can be 
referred to as the "Doubling Hypothesis" (cf. |37|, |38|). For the model at hand, this 
says that the eigenvalue expression can be written in the form 

3 

A(u, w) = Yl ai(u)$i(u, w)Ai(u), (2.9) 

i=l 

where only the "boundary contributions" (Xi(u) vary with the chosen (diagonal) K- 
matrices. The functions $>i(u, u)) and Ai(u) are doubled versions of the corresponding 
functions which appear in the eigenvalue expression for the transfer matrix with periodic 

(2) 

boundary conditions. For the A 2 model they are given by 

N 

= \\_c(u + Ui)c(u-Ui), 
1=1 

N 

$ 2 (u,u>) = X\_b(u + Ui)b{u-u)i), 
i=i 

N 

$ 3 («,w) = Y[d(u + Ui)d(u-Wi), (2.10) 
where b(u),c(u),d(u) are diagonal matrix elements of R(u) and 



A ( ) TT S ^ n ' 1 ^'- M U ^ ^ ^^[^(w — Uj) + A] 
■ jJi sinh[i(u + Uj) — A] sinh[i(w — Uj) — A] ' 

A ( ) TT s ^ n ^ L ^^ M u ^ ~ ^ sinh[|(w — — 3A] 
sinh[|(u + Uj) — A] sinh[|(u — Uj) — A] 

x tt coshl±(u + Uj)]cosh[±(u-Uj)} 



j=x cosh[i(-u + Uj) — 2A] cosh[i(w — Uj) — 2A] ' 

i r cosh[l(n + i^) - 4A] cosh[i (u - uj) - AX] 

/=iC0sh[± + -2A]cosh[±(u-«,) - 2A] ' 1 ' 

The boundary contributions a{(u) are determined by the action of t(u, oj) on the ref- 
erence state being the iV-fold tensor product of the vector with 1 as its first entry 



and elsewhere, and depend on the iT-matrices K^i^u) according to |27 



K 22e 2 K+ f I2 g 2 e 2 



a 2 (u) = K^ 2 K 22 



11 I 11 c 2 - b 2 c 2 - d 2 \ J c 2 -b 2 ))' 

K+K^e 2 Kt z K 22 g 2 K+K u g 2 e 2 



c 2 -b 2 d 2 -b 2 (c 2 - b 2 )(b 2 - d 2 ) 



«.<»> - *s k + |S - ^ " J^f) ) ■ < 2 - 12 ' 
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Here e, c, etc. are matrix elements of R(u) whose dependence on u, like those of the 
matrix elements of K ± (u), have been suppressed for convenience. For O-Sp boundaries 
we find the following: 

sinh(-u — 6A) 



ai{u) 



a 3 {u) 



c(2u) 
sinh(w — 6A) sinh(-u 



c(2u) sinh(w — 4A) 
sinh(w) sinh(-u — 2A) 
c(2u) sinh(-u — 4A) 



[u 



^{u)^-{u)t{u) 



where 



u 



2 (cosh(3A - u) ± i sinh 2A) 



(2.13) 
(2.14) 
(2.15) 

(2.16) 



and ^{u) is as defined in Eq. (|2.4j ). Together with Eqs. ( [2.10|) and ([2.1 1|) these fully 



determine the eigenvalue A(u,lj) through Eq. ( |2.9| ). The corresponding Bethe ansatz 
equations are obtained by requiring A(u, u>) to be analytic in u. It should be pointed 
out that our justification for the eigenvalue expression (|2.9|) is not as strong as we would 
like, although we have no doubt about its correctness; a derivation using the algebraic 



Bethe ansatz |29| would be preferable. In fact, a derivation using the analytic Bethe 
ansatz generalizaling the study of [^] for the quantum group-invariant case would also 
go a long way towards answering this criticism. At present our result is backed up 
mainly with numerical checks (see end of Section 2). 



2.2 Vertex and loop models with open boundaries 
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Figure 3: The "diagonal-to-diagonal" transfer matrix t£>(u) which takes the 
state («!, . . . , o>n) to the state (Pi, ... , (3^)- All internal edges are summed 
over. 

When the inhomogeneities uj are chosen to alternate as Ui = (—)' l+1 u, the Sklyanin 
transfer matrix turns into the diagonal-to-diagonal transfer matrix tv(u) for a vertex 
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model on the square lattice with open boundaries |27j (see also |31| ). More specifically, 
define the functions 



Then the object 



cd 



r b a (u) = K- b (u). (2.17) 



t D (u) = — V- (2.18) 

is the diagonal-to-diagonal transfer matrix (see Fig. 3) for the vertex model with bulk 
weights given by 




and boundary weights given by 



fa/2) \ fa/2) 



For alternating inhomogeneities u>i = {—) l+1 u we note that $2(w, w) and $3(m, cj) vanish 
due to 6(0) = d(0) = 0. Using the identity a\{u) = l\{u)r\{u) / 'R\\(u) we arrive at the 
eigenvalue expression 

t B {u)=A 1 {%), (2.19) 

with the function Ai(u) defined in ( |2.11| ). The associated Bethe ansatz equations are 
those which render A(u, uj) analytic in u, but with the inhomogeneities subsequently 
specialized to Uj = (— )i +1 u. Note that the order is crucial; analyticity of t^(u) does 
not give the Bethe ansatz equations. With this we accomplish step (I). 

We shall presently make a change of variables u — > 2iu, A — > iA + ivr/2, which takes 
us into the critical regime of the model for real u, A. The eigenvalue expression for tj)(u) 
becomes (after a shift Uj — > Uj + in) 

" smh[|(m + Uj ) + iA] smh[|(m - Uj ) + iA] 
DW sinh[i(m + Wi ) - iA] sinh[i(m - Uj ) - iA] 1 ' ' 
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while the Bethe ansatz equations now take the form (k = 1, . . . , M) 

sinh(iA — Uk) ( sinh[|('U& — 2iA — iu)} sinh[^(-Ufc — 2iA + iu)] 
sinh(iA + Uk) \sinh[^(-Ufc + 2iA + iu)} smlif^ (w^ + 2iA — iu)} 



N 



sinh[|(M fe + uj 



2iA] sinh[i(w fc — Uj) — 2iA] 



l k sinh[i(u fc + Uj) + 2iA] sinh[i(u fc - Uj) + 2iA] 
+ iA] sinh[^(-u fc — Uj) + iA] 



sinh[^(-u fe + Uj) 
sinh[i(w fc + Uj) 



iA] sinh[i(-u fc 



iA]" 



(2.21) 
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Figure 4: Allowed vertices for the loop model with partition function Z\ oop . 
Vertex i carries Boltzmann weight pi. 



Using an argument very similar to that for 0-0 and Sp-Sp [E7[] one can show that 
the vertex model partition function is equivalent to that of a loop model on the same 
lattice 



v loop 



Erf 

Q 



(2.22) 



where the sum is over all configurations Q of non-intersecting closed loops which cover 
some (or none) of the edges, with non-zero Boltzmann weights pi shown in Figure 4. 
In the configuration Q, Ni is the number of occurences of the weight of type i while P 
is the total number of closed loops of fugacity n. The explicit expressions for p« and n 
are given (in the shifted variables u, A) by 



pi = [sin(3A — u) sin(w) + sin(2A) sin(3A)] [sin(3A — u) sin(2A — u) 

P2 = P3 = ei sin(2A) sin(2A — u)~ , 

Pa = Pb — e 2 sin(w) sin(2A) [sin(3A — u) sin(2A — m)] 1 , 

P% = Pi = sin(tt) sin(2A — u)^ 1 , 



i-i 
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P8 = 1, 

p 9 = — sin(A — u) sin(it) [sin(3A — w) sin(2A — u)] 1 , 

Pio = eie -h \ 

Pll = e- iu sin[i(3A + M)]sin[i(3A-M)]- 1 , 

PX2 = eie iu , 

Pl3 = e iu cos[±(3A + u)]cos[f(3A-u)]~\ 

n = -2cos(4A) (2.23) 

where ef = e| = 1. This concludes step (II). 

To accomplish step (III) we now choose e\ = e 2 = 1 and take the "honeycomb 
limit" u — A. As seen from Eqs. ( |2.23| ), the Boltzmann weight pg vanishes in this 



limit allowing the remaining bulk vertices to be "pulled apart horizontally" to give the 
vertices of Figure 5. It is apparent that these latter vertices define a honeycomb lattice 
loop model. In other words, in the honeycomb limit the square lattice loop model 
defined through the weights of Figure 4 is equivalent to the honeycomb lattice loop 
model defined through the weights of Figure 5. There is in fact a second honeycomb 
limit, which we will consider no further, where ps vanishes and the remaining bulk 
vertices are pulled apart vertically. 

The left and right boundary weights can be normalized separately such that p 10 = 
P12 = 1, while keeping Z\ oop unchanged. In the honeycomb limit the bulk weights are 
given by p 4 = p 5 = p 6 = p 7 = p 8 = 1, p 2 = p-s = 2cos(A) = t and p\ = t 2 . Define the 
weights 

t = t = 2 cos(A), t Sp = cos(2A) cos(A)" 1 . (2.24) 

From Eqs. (|2.23|) we see that the remaining boundary weights are pn = to and P13 = tsp- 



An inspection of Figure 5 shows that we can achieve the same effect of the weight 
assignment pi by letting a honeycomb vertex (bulk or boundary) through which a bond 
runs carry a weight 1, an "empty" bulk vertex a weight t, an empty left boundary 
vertex a weight to an d an empty right boundary vertex a weight t$ p . This justifies the 



weight assignment ( |1.14j ) in the context of Figure 2, after a change of variables from A 
to g related by 

7r# = 27r-4A. (2.25) 

In the honeycomb limit, the diagonal-to-diagonal transfer matrix t-£,(u) of Figure 3 
becomes the transfer matrix % fc> r taking a row of edges up the honeycomb strip of 
Figure 1 to the next similar row. The eigenvalue expression ( |2.19| ) then becomes (with 
a rescaling Uk — > 2uk) 

A R = f\ sinh ( M i + l iA ) sinh(^ - |iA) 
M sinh(-Uj + |iA) sinh(Mj — |iA) 
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Figure 5: The vertices of Figure 4 after the honeycomb limit has been taken. 



The associated Bethe ansatz equations become (k = 1, . . . , M) 







O-Sp 



{u k , A) 



sinh(w fc — |iA) sinh(-u fc — |iA) 
sinh(w fc + iiA) sinh(-u fc + |iA) 



l N 



sinh(-Ufc + uj — 2iA) sinh(-Ufc — Uj — 2iA) 



M 

R. sinh(-Ufe + -Uj + 2iA) sinh(-Ufc — itj — 2iA) 
sinh(-u fc + Uj + iA) sinh(-u fc — Uj + iA) 

X 



sinh(uA; + uj — iA) smh(uk — Uj — iA) ' 
with the "boundary function" 9o-sp(w, A) defined as 

sinh(iA — 2u) 



6 



O-Sp 



u,X) 



(2.27) 



(2.28) 



sinh(iA + 2u) 

In the language of Figure 2, the boundary weights for the O-Sp model are given by 

ti = to, t r = ts p . (2.29) 



The corresponding weights for the 0-0 model [19], |27 

U = to, 



arc 



t Q . 



(2.30) 



The transfer matrix eigenvalue expression Q2.26| ) remains valid while the Bethe ansatz 
equations (|2.27|) have the boundary function 6o-s p (m, A) replaced by 



Qo-o(u, A) 



sinh(w — iiA) 
sinh(-u + ^iA) 



(2.31) 
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For the Sp-Sp model 



| 27j the corresponding weights are 



t,=t : 



Sp, 



t r — t 



Sp • 



(2.32) 



again with transfer matrix eigenvalue expression ( |2.26| ) and Bethe ansatz equations 
(|2.27|) with boundary function 



6 



Sp-Sp 



>,A) 



in place of 6o-s P (w, A). 
Note the relation 



cosh(u — T^iA) 
cosh(w + iiA) 



(2.33) 



6 



O-Sp 



(u, A) 2 = Q -o(u, A)e Sp -sp(^, A) 



(2.34) 



between the various boundary functions. The choice of boundary weights - (|2.29|) , (|2.30|) 
or Q2.32|) - and corresponding boundary function - Q2.28|) , (|2.31| ) or (|2.33| ) - which arise 
from the appropriate combination of A^-matrices governs the surface critical behaviour 
of the model. In the next section both universal and non-universal quantities will be 
derived from the eigenvalue expression ( |2.26|) and Bethe ansatz equations ( |2.27|) and 
their analogues for 0-0 and Sp-Sp boundaries. 

We first solved the Bethe ansatz equations numerically and compared the eigenvalues 
so obtained to those from an exact numerical diagonalization of the transfer matrix for 
the vertex model (and also the loop model). The objective was twofold; firstly to confirm 
that the results obtained using the Doubling Hypothesis ( |2.9| ) are correct, and secondly 
to obtain the Bethe ansatz root distributions corresponding to various eigenstates of 
interest. 

Restricting ourselves to even lattice sizes N, our chief findings can be summarised 
thus: Let i = N — M label the sectors of the (vertex model) transfer matrix and 
be the leading eigenvalue in that sector. Then I takes on values 0, 1, . . . , N — 1. For 
instance, A is the largest eigenvalue in the ground state sector. For all cases except 
£ = for Sp-Sp boundaries, the eigenvalue A £ corresponds to allN—£ Bethe ansatz roots 
being real and strictly positive. For the dominant eigenvalue Aq P_Sp there are N — 1 
positive roots and a pure imaginary root at approximately u± ~ i (j + £j . The scaling 
dimensions associated with are argued to give rise to the exponents characterizing 
the critical behaviour of the watermelon correlators. 
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3 Analysis of the Bethe ansatz equations 



At criticality, the free energy per site of the model is expected to scale with lattice size 



N as (see, e.g. Ref. |2T] and references therein) 



f(0 fM r 
/» = /» + ^ + ^-^ + 0(iV-3), (3.1) 

where c is the central charge for the underlying conformal field theory and here ( is a 
lattice-dependent scale factor. The bulk free energy f^, the surface energies fg and 
associated with the left and right boundaries, and the central charge c are all calculable 
from the Bethe ansatz equations, using a technique which is by now fairly standard (cf. 
En, 041 [LOL O, [12], M). In particular, the calculation for 0-0 boundaries sketched in 



19| is readily adaptable to handle the other two types of integrable boundaries. 

We are mainly interested in the eigenvalues Ah corresponding to purely real Bethe 
ansatz roots. The eigenvalue expression ( |2.27| ) and Bethe ansatz equations ( |3.4| ) are 
invariant under A — > — A and A — » A + tt. Hence it is sufficient to consider A in the 
region [0, 7r/2]. We restrict ourselves to values of A lying in the interval (0,7r/2) to 
simplify certain Fourier integrals; the expressions for the central charges and scaling 
dimensions are nevertheless valid for the end-points, as can be shown using a limiting 
procedure. The eigenvalue Ao P_Sp will be considered separately towards the end of the 
paper. Define the function 

4>(u, A) = 2 tan~ (cot A tanhw) , (3-2) 

whose derivative with respect to u is 

(f)'( A) ^ ^ 1 s * nn ( M — ^) 2sin2A (3 3) 

i du sinh(w + iA) cosh 2u — cos 2A 

After taking logarithms, the Bethe ansatz equations Q2.27D (and their analogues for the 
other boundary conditions) can be rewritten as 



N [<f)(u k , |A) + (j)(u k , | A) + Abound (wfc, A) = 2nl k + 



M 

[4>( u k + Uj, 2A) + <f)(u k - Uj, 2 A) - <p(u k + uj, A) - <p{u k - u h A)] . (3.4) 

Here the set of integers l k for k = 1, . . . , M uniquely specifies an eigenstate of the 
transfer matrix. The term Abound (u, A) depends on the boundary conditions chosen and 
is explicitly given by 

2(j)(u, |A) 0-0 
:,i("-A) { 4>(2u, A) O-Sp . (3.5) 
20(w + if,iA) Sp-Sp 
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We determined the values of Ik for the eigenvalues of interest numerically, i.e. 
by substituting the values of the Bethe ansatz roots Uj, determined from the numerical 
solution of the Bethe ansatz equations (|2.27 ), into the logarithmic form ( |3.4| ). We 



A° 


-0 


<— > 


h 




Sp 


<— > 


h 


A s P - 


Sp 


<— > 


h: 



expect that they can also be obtained by a careful consideration of the branch cuts 
involved in passing to the logarithmic form (see, e.g. ||). The following identification 
is found for those eigenstates corresponding to real Bethe ansatz roots: 

k, k = l,...,M = N-£, 
k, k = 1,...,M = N-£, 

k + 1, k = l,...,M = N-l (£^0). (3.6) 

For the eigenvalue Aq P_Sp we find that Ij = j for the real roots j = 2, . . . , N. 
Define the "counting function" z^iu) by 

1 M 1 
2ttz n (u) = 0(u,iA)+0( M ,|A)-— [<t>(u-u j ,2\)-(f)(u-u j ,\)) + —<$>(u), (3.7) 

j=-M 

where 

&(v) = ^boundK A) + <j>{u, 2A) - <f)(u, A) + <f){2u, 2A) - 0(2«, A) (3.8) 

contains the boundary contribution. For convenience, we have symmetrized zn{u) 
through the identification = — for k > and Uq = 0. Using the Bethe ansatz 
equations ( |3.4j ) it can be verified that for Bethe ansatz roots Uk we have 

z N {uk) = ^, (3.9) 

a property for which zn{u) is designed. Therefore, with the values of Ik given in (|3.6|) , 

_ dz N (u) 

Pn[u) = — (3.10) 

au 

becomes a root density for large N. Using this we can approximate for large A^ (and 
therefore also M) 

I M roc 

jj E /(«*)« J dv Poo {v)f{v). (3.11) 



-M 



3.1 Bulk free energy 

In the limit A^ — > oo Eq. ( |3.7| ) gives rise to the linear integral equation 

POO 

2n Poo (u) = 4>'(u, |A) + 0'(«, |A) - / dvpoo(v) - v, 2A) - - v, A)] (3.12) 
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for the root density Poo(u). This equation is readily solved by Fourier transforms. We 
use henceforth the notation 

a(x) = — I due lxu a(u) 



2tt 



oo 



oo 



a(u) = / dxe~ lxu a(x). (3.13) 

Using the result 

(»<*<!)• w 

which can be established by contour integration, we find that the root density is given 
by 

2cosh±Ax 

27rp 00 (x) = — . (3.15) 

2 cosh Ax — 1 



The eigenvalue expression ( 2.26|) can be rewritten in the symmetrized form 



sinh(AiA) £L sinh(w 1 - — |iA) 

Ah= . r )l' II • A FTT - 3.16 
smh(~iA) j=_ M smh(Mj — ±iA) 

Let us denote for convenience 

, . , , sinh(u — fiA) 

V(« =l°g— T7 FTT- 3.17 

smn(« — |iA) 

Then the free energy per site, defined as /at = — logAn/iV, becomes in the large iV 
limit 

POO 

foe = - dvp^v)^) 



oo 



sinh[(| — A)x] sinh(Aa;) 



dx — 7—- — — — — ^ — J —. (3.18) 

xsinh ±7ix (2 cosh Xx — 1) 



The last equality is obtained using ( 3.15 ) and the Fourier transform 



<>, . e^- 2A ^/ 2 sinhiAx / , , \ . , 

V>(* = ■ .1 2 (0 < A < f ) , (3.19) 

a; smh ^7ra; v y 

which is also readily established by a contour integration involving its derivative. Not 
surprisingly, the boundary terms have not played a part up to this point. Indeed, the 
bulk free energy is the same if one starts from periodic boundary conditions §| . In the 
next few sections we deal with the finite-size corrections in which the boundary terms 
are important. 
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3.2 Finite-size corrections 

Let us define now the functions 



P N (u) = p N (u) - poo(w), 



1 



M 



S N (u) = — S(u — Uk) - Pn(u), 

^ k=-M 

K{u) = 4>'(u,2\)-(f)'(u,\). 



(3.20) 
(3.21) 
(3.22) 



On differentiating (|3.7] ) and using (|3.12|) we find that Pn(u) satisfies the integral equa- 
tion 

1 r°° 

Pn(u) + tt / dvP N {v)K{u-v) 

Z7T J-oo 

1 r°° 1 
= j ^ dvK(u - v)S N (v) + — $'(u). (3.23) 

Taking the Fourier transform of both sides, gathering terms and then taking the inverse 
Fourier transform results in the following result for P^(u): 



j roc roo 

Pn(u) = - J7TY2 / dw / d vK ( v ~ w)S N (w)K 1 (u - v) 

(^Z7Tj J — oo J —oo 



(2vr) 2 iV 



dv&(v)Ki(u - v), 



(3.24) 



where Ki{u) is defined such that K\{x) — (l + i^(x)) . This can be given a different 
representation in terms of a function K 2 (u) defined such that K 2 (x) = —K(x)Ki(x): 

j POO I /"00 

^(«) = — / dwS N (w)K 2 (u -w) + — — / dv&tyK^u - v). (3.25) 
The Fourier transforms K(x),Ki(x), and thus K 2 (x) follow from (|3.14j ): 



K{x) 
K x (x) 



2 cosh 



\{n — 3A)x sinh |a;A 



sinh ^7nr 
sinh \tix 



(2 cosh Aa; — 1) sinh[(| — X)x] 



(3.26) 



Analogous to Eq. ( [3.23j ) we find the following expression for the finite-size corrections 
to the free energy from the definitions of and foo 



/jv — fo 



2N 



log 



1 — cos A 
1 — cos3A 



dvP N (v)ijj(v) 



dvS N (v)ijj(v). 



(3.27) 
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Substituting in the result (|3.25| ) for P/v(w), we find that 



fN-foo = --^-logf ^ C ° S n \ )--±- T du I dv&ivtyWKtiu-v) 



2N & \l-cos3A/ 4ir 2 N 



oo 



dvS N (v)e 2 (v), (3.28) 

where e2(u) is defined such that 

_ . , sinh iAx 

^ = ,(2coshL- 1) - < 3 ' 29 ) 
In fact, by contour integration, this can be inverted to give (A > 0) 

, . 1, /2cosh^ + v / 3\ 

^^ UoBhg-Vs J- (3 ' 30) 

Eqs. ( p.25| ) and ( |3.28| ) are key equations^ in what follows and lead directly to the various 
terms in the free energy relation ( |3.1| ). Observe that in deriving these equations we have 
not had to specify in which sector of the transfer matrix we work. If we choose to be 
in sector i then is to be interpreted as the lowest free energy in that sector. The 
scaling behaviour ( |3.1|) in sector leads to the central charge while that in sector i 
leads to the relevant scaling dimensions. 

3.3 Surface free energies 



The 0(1/ N) terms in Eq. ( |3.28| ) give rise to the surface energy Using the Fourier 



transform Q3.14 ) we can evaluate <3>'(x) for each type of boundary condition under 



consideration using the definitions ( |3.8|) and ( |3.5|) . For 0-0 boundaries we recover a 
result of [0 

= -II oJl^) 
2 VI — cos3A / 



sinh ^ cosh ^ cosh iZEz^ sinh (ze^> f 2 cogh ^ 



Ax 



- 4 / da; — —, , , ' '-■ (3.31) 

xsinh ^ (2 cosh Ax — 1) 



For O-Sp boundaries we find 



/s°- Sp = --log ( 1 ~ C ° SA X 
Js 2 6 ll-cos3A 




sinh ^ cosh (7r ~ 2A)z sinh (7r ~ 6A)x 
dx 2 , . 4 - Smh 4 . (3.32) 
x sinh y (2 cosh Ax — 1) 



^There are misprints in the corresponding equations given in for 0-0 boundaries. 
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For Sp-Sp boundaries we use the Fourier transform 

1 f 00 , ._ \x sinh \\x 

2' 



2tt 



+ if, (0<A<|), (3.33) 



which can also be derived using contour integration, to obtain 



2 y 1 — cos 3A J 

roc sinh 4? sinh ^ cosh ( ^~f )x cosh f 1 + 2 cosh 

+ 4 / dx - — ^ 1 * } (3.34) 

j-oo x smh ^ (2 cosh Xx — 1) 

Non-trivial numerical checks can be performed at this stage. From numerical transfer 
matrix calculations we find that the largest eigenvalue A of the transfer matrix at 
A = 7r/8 (i.e. the dilute polymer point n = 0) is given by 

1 0-0 

A = (2 + y/2) N \ (l + v^)" 1 O-Sp , (3.35) 
l + v^) Sp-Sp 



without further corrections in 1/N. Therefore the surface free energies at A = ir/8 are 
given respectively by 0, log(l + y/2) and 2 log(l + y/2). This is confirmed by numerical 
integration of Eqs. ( ggTj) , flSTffi ) and flOp . 



It is not hard to show that the surface energies for general A are related by 

/s°" Sp = \ (ft° + /s P ~ Sp ) , (3-36) 
which is, of course, a consequence of the relation (2.34) amongst the corresponding 



boundary functions. A comparison with Eq. (|3.1| ) suggests sharing out the surface free 
energies in the following way: 

AO) _ l fO-0 
Js — 2JS 

/s (s) = i /s sp-s p _ {337) 
3.4 Wiener-Hopf integral equation 

We move on now to the next leading order in 1/N in the finite-size correction J'n — foo- 
The Euler-Maclaurin formula gives rise to 

dug(u)S]yr(u) ~ — / du + du) g(u)pN(u) + 

d \J— oo J A J 

h (9(a) + 9( - A)) + mhfi (9 ' (A) + 9 ' ( - A)) ' (3 - 38) 
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for arbitrary g(u) analytic in [—A, A] where A is the largest root. In other words, A 
is such that zj^(A) = I ma , x /N where, on using ( |3.6|) , J max = N — £ in sector £ for 0-0 
and O-Sp boundaries, and J max = N — £ + 1 in sector £ ^ for Sp-Sp boundaries. The 
condition on A leads to the sum rulePI 



3& 



dup N (u) = zn(oo) — Zjv(A) 



A 

0\ / fl/f \ \ As. Jr^ \\ 1 

max 



2A /M A A Abound (go, A) _ i. 



(3.39) 



where the last equality follows from the definition (|3.7|) and the result 4>(oo, A) = n — 2A. 
Applying ( |3.38| ) to the first integral of Eq. ( |3.25|) we obtain 

I rco j poo 

Pn{u) - Pooiu) « — — : / &vKx{u - v)<& {v) - — I dvp N (v)K 2 (u - v) + 

47T iV J-oo Z7T J A 

L :AT 2 (u - A) - ^ — ttt ^^ - A), (3.40) 



4ttA^ zv y 247riV 2 pjv(A)' 

which is valid for it > A. Under this condition, certain terms implied by (|3.38| ) are neg- 
ligible because K 2 (u) decreases exponentially with increasing u and have been dropped. 
Applying ( |3.38|) to the last term of Eq. ( |3.28| ) and using the fact that e 2 {u) is an even 



function we obtain also 

dvS N (v)e 2 (v) « 2 jT dup N (u)e 2 (u) - ±e 2 (A) - 6 ^ {A) • (3-41) 

This term describes the 0(1/N 2 ) correction to the free energy and in the thermodynamic 
limit is governed by the behaviour of Pn{A) for large A. This in turn can be found by 
solving Eq. fl3.40| ) subject to the sum rule ( |3.39| ), as we now demonstrate. 
Define the functions 

X(u) = p N (u + A), 

f(u) = poo(w + A), 
1 r°° 

— - / dvK l {u + A-v)&{v). (3.42) 

T Z N J-oo 

The integral equation (|3.40| ) can then be rewritten as 



€[U) = 4**N 



1 r°° 1 K'(n) 

X (u) « /(„) + e(u) - - 1 Avx{ »)Mu - „) + —KM - 5^44,. (3.43) 

This can be recognised as a Wiener- Hopf type integral equation for xi u ) which can be 
solved by standard means. 



J This corrects a misprint in the sum rule for 0-0 boundaries given in uM 
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First split x(x) — X+( x ) + X-( x ) components analytic in the upper and lower 
half-planes H±, with 

X + {x) = ^- rdue ixu X (u). (3.44) 



(3.45) 



2tt Jo 

Then the Fourier transform of the integral equation fl3.43| ) can be written as 



X(x) « f(x) + e(x) - x+(x)K 2 (x) + —K 2 (x)C(x), 

lix 



where we define 



C{x) = -^-4 



(3.46) 



2N 12N 2 p N (A)' 

Let F(x) = f(x) +e(x) be split into components analogous to x(x). Then ( [3.45| ) can 
be rewritten as 

X-{x) + [l + K 2 (x)] [x+(x) - ±C(x)] « F + (x) + F.(x) - ±C(x). (3.47) 
Now, if we can factorize 



1 + K 2 (x] 



-i 



G+(x)G_(x) 



(3.48) 



where G + (x) (G_(x)) is analytic in Il + (respectively, n_), and G^(x)F + (x) can be split 
into + and — components 



GJx)F + {x) = Q + (x) + QJx), 
then Eq. ( |3.47| ) can be written in the form 

x+i x ) ~ 



(3.49) 



GJx) 



Q+(x) = Q-.(x)-G4x)[x-(x) + ±C{x)-F-{x) 
= P(x). 



(3.50) 



Since the first two terms are each analytic in Il + and Il_ respectively, P(x) must be a 
regular function and can be determined, for instance, from the large x behaviour of the 
first term. 

Let us now confirm the decompositions ( |3.48|) and (|3.49|) . Firstly, we find from 
( |3T26D that 

-i (2 cosh Ax - l)sinh[(vr - 2X)x/2] 



1 + K 2 (x) 



sinh(7rx/2) 



(3.51) 



It indeed factorizes as ( |3.48| ) with |T(J |TT 

GJx) -- 



2^7r(7r-2A)r (l - if) e h ^ 



r(i-f( 7 r-2A))r(i-i|x)r(|-i|x 

G-(-x), 



(3.52) 
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where 



h(x) = — 



1 71 
log- 



(7T-2A) 



log 



7T-2A 

A 



As \x\ — > oo in Il + we find, on using Stirling's formula, that 

G+(x) - 

where the coefficients are given by 



l + ^ + ^| + 0(x- 3 ) 

X X A 



(3.53) 



(3.54) 



fji 

92 



1 
3 

hi 



71 



71 



2 2(vr - 2A) 6A 



Next, from the definition ( [3.42 ) and the result (|3.15|) we have 



-ixA 



2cosh(Ax/2) 



2tt (2coshAa;-l) 



(3.55) 



(3.56) 



whose leading pole in Il_ is at x = — i7r/3A. Keeping only this leading pole term we 
obtain 



/+(*) 



e- i:i ' A 2v / 3cosh(Ax/2) 



F,(x). 



(3.57) 



2n (n — 3i\x) 

Here we have ignored the contribution to F(x) from e(x) (cf. Eq. (|3.42| )), which con- 
tributes to the finite-size corrections at a higher order in 1/N. One can now see that 
G-(x)F + (x) has only one pole in at x = — i7r/3A and use this fact to obtain 



Q+(x) 



27r(7r — 3iAx) 



(3.58) 



In particular, Q+(x) — > as \x\ — > oo in In the same limit, X+i x ) from the 
requirement that x+( u ) is integrable at the origin. Therefore from Eqs. fl3.50|) , ( |3.46| ) 
and (|3.54 ) we obtain 

Knowing C(x), G + (x), Q+(x) and P(x), we are then in a position to obtain x+( x ) 
using Eq. ( |3.50| ). We now note from the definitions ( |3.42j ) and ( |3.44j ) that x+(x) = 
/ °° due lxu Pn(u + A). Therefore from Eq. ( |3.50| ) we obtain 



-ixA 



2tt 



due ixu p N {u) = G+(x) [P(x) + Q+(x)} + —C{x). 

171 



(3.60) 



From this equation, two key results (|3.64|) and (|3.67|) will be obtained. We note at 
this point that since we have dropped e(x) - which contains the contribution from 
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$(it) and therefore Abound [u) - from F(x), Eq. ( |3.60|) does not depend on the boundary 
conditions. 

Taking the inverse Fourier transform of x+i x )i we obtain pv(A) as a contour integral, 
which can be written residue at infinity: 



x 



xG+(x) [P( x )+Q + ( x )] + — C(x) 



(3.61) 



We now subsititute Eqs. (|3~46j ), ( gjgg) , ( ^54|) and (|3T58|) into (PlflD and evaluate the 
limit to find 



7T 

3A 



(3.62) 



2 l2Np N (A)_ 

In arriving at the final result we have eliminated G + (i7r/3A) in favour of Q+{0) using 
( 13381) . Define 

p = Np N (A), 



a = N[Q + (0)+P(0)]. 
Then Eq. ( |3.62|) can be succintly written as a quadratic equation for p 



P 



m , 2tt 2 



7T 



" T + 3A Q+ 6A + 



#1 - #2 lTrpi 



12 



36A 



(3.63) 



(3.64) 



We will see in the next section that the 0(N~ 2 ) correction to the free energy can be 
expressed in terms of p, a and g\. We will also see that on using the quadratic equation 
( |3.64| ) this expression becomes just a function of a 2 . 

Let us now relate the evaluation of a to the sum rule ( |3.39 ). On setting x = in 
Eq. (|3.60|) we obtain 



dup N (u) = 2vrG + (0) [P(0) + Q+(0)] + C(0). 



(3.65) 



From ( p.52| ) we have 



G+(0) 



2J7T(7T-2A) Jh(u-2\) 



r i r 



7T 



(3.66) 



and from ( |3.46| ) we can evaluate C(0). A comparison of fl3.65| ) with the sum rule (|3.39| ) 
then leads to the result 



Aix z a 



2^,2 



2A 



7T 



-1 r 



iV-^ + -(M-iV)--+ 0bouild o (OO ' A) 
2 7r n 2rr 



(3.67) 



where the boundary contribution has finally reappeared. 
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3.5 Central charge and scaling dimensions 



We are now ready to return to the finite-size correction to the free energy ( |3.41| ) . From 
Eq. (|3.30|) we have 

e 2 (u) ~ V3e 3A 



(3.68) 



for large u. Substitution into Eq. ([3 .41]) yields 

dvS N (v)e 2 (v) w 

2\/3 [°° due m (£) p N (u) 
Ja 



V3 



V3 

N 18XN 2 p N (A) 



7rA 

e . 



(3.69) 



Using the equation ( p.60|) with x = i7r/3A and the expression Q3.46Q for C(x) we end 
up with 



dv S N (v)e 2 (v ) 4v / 3vre 3aG h 



P 



3A 



3A 



(3.70) 



Both G+ (§) and Q + (§) can be eliminated in favour of Q+(0) using ( |3.58| ). This in 
turn can be expressed in terms of a and p. We have thus 



dvS N (v)e 2 (v) 



4^ 



7T 



AT2 



7T 



2i 2 

— Q+ 3'- 



m 

2 n P/ 
i \ 



1 i<7i 
2 ~ 8^ + 487rp + 72Xp 



(3.71) 



Expanding out the terms and using the quadratic equation ( [3.64 ) for p we find eventu- 
ally that 

7T / 2 N 



dvS N (v)e 2 (v) 



24N 2 1^3, 



1 - 48vr V 



(3.72) 



Now from Eq. ( 3.67|) we have the result 



47Ta 



2„ 2 



-1 



2A _i_ 1 

7T "T 2 



0-0 

£(l-2£)_2A] Z O-Sp 

^(l-v)-v-^] 2 Sp-Sp (1^0) 



(3.73) 



for the various eigenstates of interest, except for the ground state with Sp-Sp boundary 
conditions. Let us first consider the ground states (£ = 0) for the 0-0 and O-Sp cases. 
Substituting the result for a 2 in Eq. ( |3.73| ) into Eq. ( |3.72j ) and comparing with the 
result (|3.1j) leads to the identification ( = 2/y/3 for the scale factor|, together with the 
§This corrects misprints in p9, p§| 
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central charges 



o-o 



c 



„Q-Sp 



1 - 



3(tt-4A) 2 

7t(7T - 2A) 

3(4A) 2 



7r(7T - 2A)' 

With the change of variables from A to g related by (|2.25 ), we obtain the results fll 
and Ql.15 ). The mixed-boundary scaling index t° _Sp defined by 



(3.74) 
(3.75) 



,o-s P _ c o-o _ 24 ^o-sp 



(3.76) 



is therefore given by t " Sp = (3 - 2g)/4g = h h2 in the Kac formula ( |i~T0| ). Using the 
notation q for the effective central charge for the eigenstate labelled by £, we find from 
(|3772D and (|3~73p that 



x>-o 



c 



o-o 



O-Sp O-Sp 



-24 
-24 



The scaling dimensions defined through the inverse correlation lengths [32 



-i 



log-:- ~ 



N 



(3.77) 
(3.78) 

(3.79) 



are then given by ( |1.9|) and ( |1.16|) . 

Let us now turn to the ground state for Sp-Sp boundaries, where there is an imagi- 
nary Bethe ansatz root and the preceding analysis has to be slightly modified. Firstly we 
noted numerically that the position of the imaginary root reaches its asymptotic value 
much faster than 1/N. We therefore consider it to be fixed; and that this assumption 
does not affect the finite-size corrections to the order in which we are interested. For 
convenience we perform a shift Uj — > Uj-i in the eigenvalue expression ( 2.26|) and as- 
sociated Bethe ansatz equations. The imaginary root is now uq ~ i (j + | J while the 
eigenvalue expression becomes 



A Sp-Sp 



fiA) 



sinh(u + |iA) sinh(ti 

sinh(tt + |iA) sinh(wo — §iA) M smh(uj + ^iA) smh(uj — |iA) 



|iA) ^j-j- 1 sinh(uj + |iA) sinh(wj 



(3.80) 



The Bethe ansatz equations ( |3.4| ) are now valid for k = 0, 1, . . . , M = N — 1, with 
Ik = k + 1. The product term in (|3.80|) , which involves only real roots Uj can now be 
treated as before. The only difference is that the counting function zn{u) has to be 
modified in order that z^{uk) = Ik/N holds for the real roots U}~ in the shifted variables. 
We find that <&(w) in ( |3.7|) has to be replaced by 



$§ p_Sp (u) 



+ 



0bound(w, A) - (j)(u + M , 2A) + 

(f)(2u,2X)-(f)(2u,X), 



(u + u , A) 



(3.81) 
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where unlike previously uq is no longer identified with 0. With these changes, the anal- 
ysis of finite-size corrections then proceeds as before. Since we have already obtained 
the surface energy from the excitations, we restrict ourselves to the 0(N~ 2 ) terms. We 
find that Eq. (|3.72|) is unchanged. The sum rule ( |3.39| ) on the other hand is now 



dup N (u) 



1 2A 



due to $ 



Sp-Sp, 



u 



i. This results in the replacement of ( |3.73| ) with 



4ttV 



1_ 2A 

Equations d3~72|) and ( ggg ) then give rise to 



1 _ 2A 

2 7T 



c 



Sp-Sp 



3(tt-4A) 2 
7T(7T - 2A) ' 



(3.82) 



(3.83) 



(3.84) 



which is the same as for 0-0 boundaries. From Eqs. ( |3.72| ) and ( |3.73| ) we can now 
calculate the effective central charge for A^ p_Sp : 



c s P -Sp _ c sp-s P = _ 24 



f (i-f) _1 +H-f)^+(i-^ 



(3.85) 



The scaling dimensions X^ p Sp are then given by ( |1.13 ). 



4 Discussion 

In this paper we have studied the 0(n) loop model on the honeycomb lattice with 
open boundaries from the point of view of the underlying Izergin-Korepin i?-matrix 
and its associated i^-matrices. Three inequivalent sets of integrable boundary weights 
and their corresponding Bethe ansatz solutions were thus obtained. Two of these sets 
were of "pure-type" |27J and one of mixed-type. One of the pure-type boundary weights 
(denoted here by O-O), with the bulk and surface couplings being the same, was pre- 
viously studied using the coordinate Bethe ansatz |B| and argued to correspond to 
the 0(n) model at the ordinary surface transition. By analysing the associated Bethe 
ansatz equations for the second set of pure boundary weights (denoted here by Sp-Sp) 
we found agreement with known results (some of which were conjectures) for the spe- 



cial surface transition |fl5| , pOf , which in the n — ► limit corresponds to the adsorption 
transition for polymers. This justifies our claim to have found the critical couplings for 
the 0(n) loop model on the honeycomb lattice at the special transition and allows us 
to "prove" the abovementioned conjectures. The set of mixed-type boundary weights 
therefore should correspond to the 0{n) model at the mixed ordinary-special surface 
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transition. Our Bethe ansatz analysis gives a central charge and mixed-boundary scal- 
ing index which are indeed in agreement with a recent result plj for this transition, 



providing further evidence for our claim. In addition our analysis has furnished the 
scaling dimensions for this transition. 

(2) 

It is interesting that the i^-matrices for the A 2 model have "accounted for" the 
ordinary and special surface transitions in the manner described in this paper. This 
leaves the interesting open problem of finding an integrable lattice 0(n) model which 



exhibits the extraordinary surface transition p6) which corresponds to having fixed 0(n) 



spins at the boundary, and for which conformal invariance predictions exist |T5|, ^T[]. It 
is also an interesting open problem to derive the physical boundary S-matrices, which 
constitute a starting point for the techniques of ||20|| , from the i^-matrices. This ought 
to be achievable using the Bethe ansatz analysis performed here, along the lines of j|3 
for the XXX model, or using the vertex operator approach p4 
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